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Abstract 

If dark matter possesses a lepton number, it is natural to expect the dark-matter annihilation 
and/or decay mainly produces the standard model leptons, while negligible amount of the an- 
tiproton is produced. To illustrate such a simple idea, we consider a scenario that a right-handed 
sneutrino dark matter decays into the standard model particles through tiny i?-parity violating 
interactions. Interestingly enough, charged leptons as well as neutrinos are directly produced, and 
they can lead to a sharp peak in the predicted positron fraction. Moreover, the decay of the right- 
handed sneutrino also generates diffuse continuum gamma rays which may account for the excess 
observed by EGRET, while the primary antiproton flux can be suppressed. Those predictions on 
the cosmic-ray fluxes of the positrons, gamma rays and antiprotons will be tested by the PAMELA 
and FGST observatories. 
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I. INTRODUCTION 



The presence of dark matter has been securely established by numerous observational 
evidences. In particular, the latest 5yr WMAP data determined the dark matter abundance 
with unprecedented precision as [l| 

tt DM h 2 = 0.1099 ±0.0062, (1) 

where Qdm is the fraction of critical density in dark matter and h is the Hubble parameter 
in units of 100 km/Mpc/sec. In the meantime it is not known yet what dark matter is 
made of despite many experimental direct /indirect searches hitherto — but there is a hope 
that the PAMELA Q and FGST (formerly GLAST) satellites may provide us with some 
important information on the nature of dark matter. 

It is normally assumed that the dark matter is charged under an exact discrete symmetry 
to ensure its stability, for instance, i?-parity in a supersymmetric theory. There might exist 
many discrete symmetries realized in our vacuum and one of them may be responsible for 
making the dark matter stable. However, if the discrete symmetry is broken, the dark matter 
will be unstable and eventually decay into the Standard Model (SM) particles. Recently 
such a decayingdark matter has attracted much attention (e.g. the gravitino with i?-parity 
violation 0, Q, Q, B] or a hidden U(l) gauge boson {9J and gaugino [loj]), since the high- 
energy cosmic rays produced by the decay of dark matter may account for the observed 



anomalous excesses in gamma rays ll|, |12( and/or positrons 13j. 



Very recently, the PAMELA data on the antiproton flux has been released [ijj], and 
it suggests that the observed antiprotons are mainly secondaries, which is consistent with 
the balloon-borne experiments at the top of atmosphere (TOA) [15(|. This will place tight 
constraints on possible dark matter candidates for explaining the anomalous excesses in 

n 

positrons and/or gamma rays [16j]. Of course, the predicted antiproton flux has a large 
uncertainty mainly due to our poor understanding of the cosmic-ray propagation inside 
our Galaxy. However, if we take the constraint on the antiproton flux seriously, and if 
we attribute the excesses in positrons and gamma rays to the dark-matter annihilation or 
decay, we are led to explore a dark matter candidate naturally satisfying the constraint on 
the antiproton flux. Our main idea is as follows. If the dark matter has a lepton number, it 
is quite natural to expect that it annihilates or decays mainly into the leptons, while only 
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negligible amount of the antiprotons are produced. The purpose of this paper is to illustrate 
this simple and naive idea by using an explicit example. 

There are several candidates for leptonic dark matter such as the sterile neutrino [l^]. 
As an example, in this paper we consider a scenario that a right-handed sneutrino vr is 
the lightest supersymmetric particle (LSP) and accounts for dark matter of the universe. 
The neutrino mass is assumed to be Dirac type in our study. If the -R-parity is an exact 



symmetry of nature, the vr dark matter is absolutely stable [18|, [19|, [20]. However, the 
.R-parity may not be an exact symmetry and is only an approximate one accompanied with 
tiny violations. We focus on the -R-parity violating bilinear term throughout this paper. In 
the absence of the -R-parity, the vr dark matter is not stable anymore, and directly decays 
into the charged leptons (r, /i and e) and neutrinos, and it can also decay into the quarks 
and the W and Z gauge bosons depending on the mixing with the Higgs bosons. As a result, 
a sharp peak is predicted in the positron fraction, which may explain the anomalous excess 
observed by High Energy Antimatter Telescope (HEAT) [131], MASS |2l| and AMS [22] 
experiments. Interestingly, the preliminary PAMELA data also exhibits such anomalous 
excess in the cosmic-ray positron fraction, although we do not try to include the preliminary 
data since it is not formally released yet. The continuum gamma rays are also produced 
mainly from the decay of pions, and those gamma rays may account for the excesses observed 



by EGRET 



111 ]. For a proper set of parameters, the production of antiprotons can be 



suppressed, which is consistent with an observational fact that the antiprotons measured by 
the balloon-borne experiments []jj] and also by the PAMELA satellite 14] are considered to 
be mainly secondaries 23|]. The suppression in the antiproton flux is particularly important 
because some decaying dark matter scenarios predict too large antiproton flux at the solar 
system 

This paper is organized as follows. In Sec. [Ill we present the effective interactions and the 
relevant decay modes used in our calculations. In Sec. IHIl we calculate the spectra of the 
positron, gamma-ray, and antiproton produced from the decay of the ur, and compare them 
with the observational data. We also discuss how the right-handed sneutrinos are produced 
in the early universe in Sec. IIVI and we give our discussions and conclusions in Sec. [V] 



3 



II. FRAMEWORK 



The non-vanishing neutrino masses have been firmly established by neutrino oscillation 



experiments (see Ref. 24j for recent review and references therein), although it is not known 
yet whether the neutrinos are Majorana or Dirac fermions. In this paper, we consider a 
Dirac-type neutrino mass in a supersymmetry theory without i?-parity conservation, and 
assume the gravity mediation. We introduce two additional interactions to the minimal 
super symmetric standard model (MSSM); one is the neutrino Yukawa coupling for the Dirac 
neutrino mass, and the other is the .R-parity violating bilinear term. The superpotential is 
therefore * x 

W = WWi + V-jNiLjHu + mH u Li, (2) 
Wmssm = y?UiQiH u - viD % QiU d - yfEiLiHi + [iH u H d , (3) 

where N is the right-handed neutrino superfield, y u is the neutrino Yukawa coupling, /ij 
denotes the coefficient of the i?-parity violating bilinear term, and the indices % and j denote 
the generations. We neglect the flavor mixings in the MSSM Yukawa interactions, and 
assume the minimal Kahler potential for all the MSSM particles as well as the right-handed 
neutrinos. 

Let us first discuss the neutrino Yukawa coupling. The three left-handed neutrinos in the 
weak eigenstate z/, are related with the mass eigenstates z), as 

Vi = Uij ■ Oj, (4) 

where i runs from 1 to 3 and the well-known mixing matrix U takes the following form: 

~~ Sl2 c 23 — S 23 Si3Ci2 C23C12 — S23S13S12 S23C13 , (5) 

\ S23S12 - S13C23C12 -S23C12 - S13S12C23 C23C13 j 
where c^- = cos 9ij and Sy = sin 0^ , and 9^ denotes the mixing angle of neutrinos z/j and Vj . 
We can define the right-handed neutrinos so that the neutrino Yukawa coupling matrix is 



C13C12 S12C13 S13 



^ 1 The right-handed sneutrinos acquire tiny vacuum expectation values (vev) after the electroweak symmetry 
breaking in the presence of the bilinear i?-parity violating term. However, the expectation values are 
suppressed by the neutrino Yukawa couplings, and so, we can neglect their effects. On the other hand, 
the bilinear term can be induced from the neutrino Yukawa couplings with non-vanishing vevs of the 
right-handed sneutrinos or the linear terms in the Kahler potential. 
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Am§! [eV 2 ] 


lAmfil [eV 2 ] 


sin 2 9 \2 


sin 2 6*23 


sin 2 6>i3 


best fit 


7.651^ x 10- 5 


2.40±g;^f x 1(T 3 


n 304+0-022 


o ™ +0 - 07 


01 + 016 
u - Ui -o.on 


adopted 


7.65 x 1(T 5 


2.4 x 10~ 3 


0.304 


0.50 






TABLE I: The best-fit and adopted values of three-flavor neutrino oscillation parameters from 
global data, including solar, atmospheric, reactor (KamLAND, CHOOZ) and accelerator (K2K) 



25|. 



experiments 
given by 

y v = diag(mi,m 2 ,m 3 )f/ t , (6) 
where m; is the mass of z)j. In Table HJ we show the observational constraints on those mixings 



and the mass differences based on the global three- neutrino analysis [25( . Since the neutrino 
oscillation data is not sensitive to the absolute masses of the neutrinos, the following three 
neutrino mass spectra are possible: (i) normal hierarchy case m 3 > m 2 ^> m\ (ii) inverted 
hierarchy case m2 > m\ > and (iii) degenerate case mi ~ irt2 — m?,. Throughout this 
paper, we adopt the normal hierarchy for simplicity, therefore, m\ ~ Am^, m| ~ 
and we adopt massless mi in our numerical study. 

Next we consider the .R-parity violation. In the presence of the bilinear .R-parity violation, 
it is known that there is no unique way to define Lj and since they have the same quantum 
numbers 26|. Taking account of the soft SUSY breaking terms, the left-handed sneutrino 
(L°) acquires a non-vanishing vacuum expection value (vev) after the electroweak breaking 
in a general basis. In our work, we adopt a basis such that the (L°) vanishes by a proper 
redefinition of Li and Hd- In this basis, the trilinear i?-parity violating interactions are 
generically induced. Furthermore, the right-handed sneutrinos get mixed with the left- 
handed sneutrinos as well as up- and down-type Higgs, which would make analysis on the 
right-handed sneutrino decay complicated. The purpose of this paper is not to explore all 
possible parameter space, but to illustrate our basic idea that the dark matter with a lepton 
number can account for the sharp rise in the positron fraction while the antiproton flux can 
be suppressed. Thus we focus on a case that the dark matter is comprised of N%, which 
decays into the SM particles through the i?-parity violating bilinear term with fii ^ 0. 
Then, the relevant mixings of with the left-handed sneutrinos and the Higgs bosons are 
suppressed by the small mixing angle #13. In particular, those mixings are absent in the 
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limit of #i3 = 0, which therefore greatly simplifies our analysis. The adopted values of the 
mass differences and the mixing angles are shown in Table [B We will discuss later how our 
result is modified for other choices of the parameters. 

In our set-up, the main decay channels of the right-handed sneutrino, z>^ 3 (the scalar 
component of iV 3 ), are neutrinos and charged leptons through the neutrino-neutralino and 
charged-lepton-chargino mixings. The corresponding Feynman diagrams are shown in Fig. [TJ 
We can see from Eq. ([3]) that the interactions between the right-handed sneutrino vr with 
the higgsino and the SM leptons are given by 



(7) 



where v c Ri denotes the scalar component of iVj, and i\ = e, £2 = \i and £3 = r. First let us 
consider the neutrino-neutralino mixing. The mass matrix of the neutral fermions Mjv is 
given by C D (-1/2)^ 0T M N ^° with [26| 



M 



M (N) M (N) 
1V1 MSSM IV1 R 



N 



M 



(N)T 



(8) 
(9) 



T , Mjrfg SM is the usual neutralino mass 



where we have defined \1/ = (B, W 3 , H%, H®, v L 
matrix in MSSM, and Ml is the mass matrix between the neutralinos and the left-handed 
neutrinos arising from the R— parity violating interactions. The explicit expressions of 



M MSSM and M R ] are 



M 



(N) 

MSSM 



Mi 
M 2 
-m z s w C(3 m z c w c p 
\ mzswsp -m z c w sp 



-m z s w cp m z s w si3 
m z c w cp -m z c w sp 
-fj, 



\ 







M 



(TV) 




/ 



(10) 

where cw = cos Qw and Sw = sin 9 W with Q w being the weak mixing angle; cp = cos (3 and 
sp = sin/3 with tan/3 = v u /vd where v u (vd) is the vev of up-type (down-type) Higgs field. 
Here we keep /i2 and /X3 in the mass matrix, although we will set them to be zero in the 
following analysis. The mass matrix Mn could be diagonalized by an unitary matrix V, 
and the mixing between the up-type higgsino H® and the left-handed neutrinos v Li is given 
by V4i+4. In a similar way, the charged leptons mix with the charginos via the i?-parity 



6 



(a) 



j 

(b) 



FIG. 1: Feynman diagrams of right-handed sneutrino decay via the mixing between the up- type 
higgsino and the SM leptons. 

violating interactions. The mass matrix Mc of the charginos and the charged leptons is 
given by C D -^~ T M C ^ + with [26] 



(c) 

Mc=| M 2x3 

Mf ] y°5 ijVd /V2 

where we have defined Vp- = (W~, H~ 4 t ) T and ^ + = (W + , H+, ef) 7 ', M$ SSM is the usual 

(C) 

chargino mass matrix in MSSM, M y R 1 is the mass matrix induced by i?-parity violating 
interactions and y e is the charged lepton Yukawa coupling. The explicit expressions are 



(O M 2 V2m z c w s p | {c) 

m MSSM ~ ' 1 E 



/ -//! ^ 



12) 



-/x 2 

-/i 3 y 

The mass matrix Mc can be diagonalized by two rotational matrices 0_ and + which 
relate the gauge eigenstates of and with their mass eigenstates, respectively, i.e. 
7 McO + = diag (m--,m--,m e ,m /J ,m r ). The mixings between H£ and the right-handed 
charged leptons are given by the elements [0 + ] 22+i . Combining the interaction (J7J) and 
the mixings between the up-type higgsino and the SM leptons, the vr will decay into the 
SM leptons as shown in Fig. [TJ 

Finally, we summarize the relevant decay processes for our study below; 

v Rl - £je c + oc \y% [0 + ] 2i2+fc | 2 , (13) 

where V and 0± are the rotation matrices in the neutralino-neutrino and the chargino- 
lepton respectively. In Table [Til we show the branching ratios of a decaying z/r 3 , in which 
we assume non-vanishing with /x 2 , 3 = for simplicity. Notice that positrons are always 



channel 


e + p 


e + r 


vv 


branching ratios (%) 


1.5 


1.5 


97 



TABLE II: Branching ratios of right-handed sneutrino v^. 

produced in the charged lepton decay modes since only \i\ is non-vanishing. Also, because 
of the smallness of the positron mass, the decay branching ratios of positron are highly 
suppressed down to a few percentages, compared to that of the neutrino production. 



III. COSMIC-RAY FLUXES 



Let us here mention that both the right-handed sneutrino and its antiparticle can be 
dark matter and contribute to the cosmic-ray signals. For simplicity, we assume that both 
of them have been produced with an equal amount in the early universe, therefore, we have 

n DM = =2^ = 2Q, R , (14) 

Pc Pc 

where Qdm and Q„ R are density parameters of the dark matter and the right-handed sneu- 
trino, respectively; p c is the critical density and po R {p c ) is the energy density of right-handed 
(anti) sneutrino. In the rest of this paper, contributions of both right-handed sneutrino and 
its antiparticle are included in the cosmic-ray fluxes. The method for calculations of the 
gamma-ray flux and positron fraction are the same as that in our previous study 9|, there- 
fore, we only show the equations which are needed in the calculations. We refer readers who 

Q 

are interested in the derivations to Ref. [7|] and references therein. 



A. Positron fraction 



After being produced from the decay of the right-handed sneutrino, positrons will prop- 
agate in the magnetic field of the Milky Way. The positron flux is given by 

^ m ^=^r["' l2iE ' G{E ^ )i ^ (i5) 

where E is in units of GeV, and m i>R and Tp R are the mass and lifetime of the right-handed 
sneutrino dark matter. Here dN e + / dE is the energy spectrum of positron from the dark 
matter decay. We have used PYTHIA [27] to calculate the energy spectra, and the numerical 



S 



0.1 



0.01 - 



LU 



0.001 



0.0001 - 



le-05„ 



1 1 II 1 II 1 1 II II II 1 II 


mm i i- 




i i i 1 1 1 n| i i i 1 1 1 in i i i 1 1 1 1 1 


1 1 


Mass - 1 00 GeV 






— e + Mass - 300 GeV 


- 






0.01 


_ V 










- V 

V 

T 








0.001 




- 






: 


- 






0.0001 






i i i ill 


1 i ml i i 


le-05 


1 1 <l I 


1 1 



0.1 



1 10 
E (GeV) 



100 



0.1 



1 10 
E (GeV) 



100 



FIG. 2: Energy spectra of gamma, positron and neutrinos generated from the decay of a vr. 
results are shown in Fig. [2j The Green function G(E, E') in Eq. (|T5!) is approximately given 

by & 



G(E,E'] 



1 ni6 

^ 1U r a+b(E s - 1 -E' s - 1 )nr m/ 



E 2 



^9{E' - E) [sec/cm 3 



(16) 



where 5 is related to the properties of the interstellar medium and can be determined mainly 
from the Boron to Carbon ratio (B/C) 28j|. In the numerical study, we adopt the following 
parameters, 5 = 0.55, a = —0.9716 and b = —10.012 [7|], that are consistent with the B/C 
value and produce the minimum flux of positrons. 

In addition to the positron flux from dark matter decay, there exist secondary positrons 
produced from interactions between the primary cosmic rays and nuclei in the interstellar 
medium. The positron flux is considered to suffer from the solar modulation, especially for 
the energy below 10 GeV. If the solar modulation effect is independent of the charge-sign, 
one can cancel the effect by taking the positron fraction, 



(17) 



which is indeed measured in many experiments. To estimate the positron fraction, it is nec- 
essary to include the electron flux. We use the approximations of the e~ and e + background 
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FIG. 3: The fraction of positron flux 
together with experimental data 



13, 



21 
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where E is in units of GeV. Therefore, the fraction of e + flux is 



prim 
e+ 



prim 



Qsec 
e 



(18) 



(19) 



where k is a free parameter which is used to fit the data when no primary source of e + flux 



exists 



30 



Our numerical results for rriy = 100 GeV and 300 GeV are shown in Fig. 



The typical feature is the sharp turnup at E « 10 GeV and a drop-off at E = due 
to the contribution of the direct production of e + from the decay of dark matter [9j . 



B. Gamma-ray flux 

Generally, the main contribution to the gamma-ray spectrum arises from the pion 7r° 
generated in the QCD hadronization process of qq and in the decay of r. Since no quark 
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pair is produced in our set-up, the only source of 7r° is from the decay of r. To estimate the 
spectrum, we use the PYTHIA 2t| Monte Carlo program with the branching ratios shown 
in Table HT1 and the energy spectra dN 1 /dE are presented in Fig. [21 It is worth noting that 
there is no line emission of the gamma rays from the decay of Ur, which is present in the 
case of the gravitino dark matter 0, Q]. 

There are galactic and extragalactic contributions from the decay of vr to the observed 
gamma-ray flux. The flux of the gamma-ray from the extragalactic origin is estimated as 

in 

_E 2 cQ DM p c 

eg 



dE 



.1/2 



y ^ , dN~ 



V 



-3/2 



Ha 



(20) 



where c is the speed of light; Qm and Q\ are the density parameters of matter (including 
both baryons and dark matter) and the cosmological constant, respectively; H is the Hubble 
parameter at the present time; y = 1 + z, where z is the redshift, andy eq denotes a value of 



y at the matter-radiation equality. For the numerical results, we use [l] 

Pt 



Vt DM h 2 = 0.1099, n M h 2 = 0.1326, fi A = 0.742, p c = 1.0537 x 10" 5 GeV/cm 3 . (21) 



On the other hand, the gamma-ray flux from the decay of dark matter in the Milky Way 
halo is 



dE 



E 2 



dN~, 



halo 



^nmy R Ty R dE 



phaMdi 



l US 



(22) 



where phaio is the density profile of dark matter in the Milky Way, (^j los phaio{$)dtj is the 
average of the integration along the line of sight (los). In our calculation, we adopt the 
Navarro- Frenk- White (NFW) halo profile [32] 

Po 



p(r) 



(23) 



(r/r c )(l + r/r c ) 2 ' 

where r is the distance from the center of Milky Way, r c = 20 kpc, and po is se t in such a way 
that the dark matter density in the solar system satisfies p(r ) = 0.30 GeV/cm 3 33] with 
r = 8.5 kpc being the distance from the Sun to the Galactic Center. For the background, 
we use a power-law form adopted in Ref. 

- bg 

where E is in units of GeV. The predicted gamma-ray flux from the decaying vr are shown 
in Fig. H] together with experimental data. Since no hadronic decay mode exists and the 
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FIG. 4: Gamma-ray flux predicted from decay of right-handed sneutrino dark matter, shown 
together with the EGRET data [11]. 



amount of r is insignificant, the predicted gamma-ray flux is below the EGRET data for the 
parameters we have chosen. However, we do not regard this tension as a serious conflict, 
and we expect that the upcoming data from FGST should be able to further probe such a 
prediction from the right-handed sneutrino dark matter decay. Also, we should emphasize 
that the relatively suppressed gamma-ray flux shown in Fig. H] is due to the simplification 
we have adopted. The gamma-ray flux may fit the EGRET data if we allow // 2 ,3 to be 
non-zero and/or the other generations of the right-handed sneutrino to contribute to the 
dark matter abundance. This is because the decay branching ratios of r will increase and 
qq decay channel will open due to the mixing between the right-handed sneutrinos and the 
Higgs bosons, then obviously, more ir° will be generated. 



C. Antiproton flux 

As we have mentioned before, the predicted antiprotons from decay of gravitino dark 
matter Q] tend to have a tension with the observational data. Of course, this apparent 
tension may not be serious at all, considering the large uncertainties on both experimental 
data and diffusion models. However, we should better have a scenario in which the antipro- 
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ton flux could be suppressed. In our set-up, the right-handed sneutrino does not decay into 
quarks (at least the decay rate is suppressed), and therefore the antiproton flux is negligible. 
This feature will be lost once we allow # 13 and/or /x 2 ,3 to deviate from zero. Or if Vr\^ a ls° 
contribute to the dark matter abundance, their decay generically produces the antiproton, 
since they can mix with the Higgs bosons as well as the left-handed sneutrinos. The antipro- 
ton flux in generic parameter space will depend on those parameters as well as the neutrino 
mass spectrum, and it is beyond the scope of this paper to survey all the possibilities. How- 
ever, it is important to keep in mind that there is a set of parameters where the antiproton 
production becomes negligible. 



D. Neutrino flux 



The main decay channel of the is vv as can be seen from Table [III It has been 
recently studied in detail whether the neutrinos produced by the gravitino decay can be 



detected by the current and future observations [34J. In our scenario, the neutrino flux is 
enhanced compared to the positron and the gamma-ray ones, and therefore the neutrino 
signal from the decay can be larger by an order of magnitude compared to the gravitino 
of the same mass. Thus, the neutrino signal, especially u T , may reach the sensitivity of future 
experiments, e.g. Hyper-Kamiokande. However, we need to mention that the enhancement 
in the neutrino production will be lost for another choice of the parameters. 



IV. COSMOLOGICAL PRODUCTION OF v R 

We would like to briefly discuss here how the right-handed sneutrino is produced in the 
early universe. In order to account for the observed dark matter abundance, the cosmological 

n 

abundance of vr must satisfy [1| 

tt^h 2 + tt^h 2 = 0.1099 ±0.0062, (25) 

or equivalently, 

n~j„ + n,~,c 

4.2 x 10 



where n i>R is the number density of z/r, and s the entropy density of the universe. If not 
only the lightest but also heavier right-handed sneutrinos are stable in a cosmological time, 
we need to sum over those right-handed sneutrinos in Eq. fl26l) . 
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In the case of the neutralino dark matter, it is usually assumed that the neutralinos are 
thermally produced at a temperature above the freeze-out temperature. However, since the 
neutrino Yukawa coupling y v is very small (~ 10~ 13 ) in our scenario, the vr is never in 
equilibrium unless it has other unknown strong interactions with the radiation in the early 
universe. Therefore we need to consider non-thermal production of vr. There are several 



possibilities [18[, and we consider one by one as follows. 

First of all, vr is produced by the decay of the lightest SUSY particle in the MSSM 
sector (MSSM-LSP) through the neutrino Yukawa coupling(s). In the presence of .R-parity 
violating interactions, the MSSM-LSP also decays into the SM particles without producing 
vr. If the .R-parity violating coupling is much larger than the neutrino Yukawa coupling, 
the production of vr will get suppressed. This is the case in our set-up, since a relatively 
large -R-parity violation is needed to compensate the suppression of the decay rate into 
the charged leptons. However, for another choice of the model parameters, the magnitude 
of the i?-parity violation tends to be as small as the neutrino Yukawa coupling, and the 
MSSM-LSP decay into vr via the -R-parity conserving interactions may occur at a non- 
negligible rate. This would be indeed the case if the neutrino masses are degenerate, since 
the neutrino Yukawa coupling will become larger than that in the case of normal hierarchy, 
and the needed coupling strength of the -R-parity violation to account for the positron excess 
becomes correspondingly smaller. Then some amount of vr will be produced by the decay 
of the MSSM-LSP. In the .R-parity conserving case, the production from the MSSM-LSP 



was studied in Ref . [18] . We can estimate the abundance in a similar way, taking account of 
the suppression due to the presence of the decay via -R-parity violating interactions. That 
is to say, contribution to the density parameter of vr from the decay of the MSSM-LSP, \i 
is given by 

AQ, R = B u Q x , (27) 

\m x J 

where B u denotes the branching ratio of the \ decay into ur, and m x is the mass of \. 
Note that, due to the presence of the -R-parity violation, the x decays faster than the case 
without it, which makes it easier to satisfy the big bang nucleosynthesis (BBN) bound on 
the MSSM-LSP decay. However, it also means that larger abundance of the MSSM-LSP 
is needed to produce a right amount of vr. In the normal hierarchy case, this would push 
up the soft masses of the SUSY particles, and the naturalness may be called into question. 
On the other hand, in the case of degenerate neutrino masses, the effect of the -R-parity 
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violation is relatively smaller, and we expect that a right amount of vr can be generated by 
the MSSM-LSP decay 

Another interesting possibility is that the gravitino decay produces vr. The gravitinos 



are generated by particle scatte ring in thermal plasma 35l . |36| | , and they are also generically 



produced by the inflaton decay 37|, |38|, |39j, |40|, |4l|. For the gravitino mass is of 0(100) GeV 



it typically decays into the SM particles and their superpartners during BBN. The energetic 
decay products can significantly change abundances of light elements, thus the gravitino 
abundance is tightly constrained by BBN. If the gravitino is the next-to-lightest SUSY par- 
ticle (NLSP), however, it will mainly decay into vr and a neutrino, and therefore, the BBN 
constraint will be greatly relaxed. To account for the abundance Eq. (1261) . the reheating 
temperature should be as high as about 10 10 GeV if the gravitino is mainly thermally pro- 
duced. Note that the production of vr from the gravitino decay can be concomitant with the 
production from the MSSM-LSP decay as well as the following mechanisms, so the reheating 
temperature can be lower. 

The inflaton decay may also directly produce vr. Indeed, if the inflaton acquires a non- 
vanishing vev at the potential minimum, it couples to all the matter fields as well as the 



gauge fields in supergravity 



39 



40 



41] . However, the smallness of the neutrino Yukawa 
coupling makes the branching ratio of the decay into vr extremely small. It is likely that we 
need to assume relatively strong interactions between vr and the inflaton sector (or other 
fields whose energy dominates the energy density of the universe after inflation) for the 
enhancement of contribution from inflaton decay. 

The last possibility is the scalar condensation of vr. If the position of vr is deviated 
from the origin during inflation * 2 , Vr will start oscillating after inflation when the Hubble 
parameter becomes comparable to its mass. The abundance is estimated to be 

^ Tr_ (vj?\ 2 / mgg W / T R \ ( ujt V . , 

s A mi/R \M P ) ~ UOOGeW Vl0 6 GeV; Vl0 10 GeV; ' 1 ' 

where Tr denotes the reheating temperature of the universe and v° R sc is the amplitude of 
the oscillations. We have assumed that the sneutrino starts to oscillate before the reheating 
is completed. For an appropriate choice of the amplitude, it can account for the right 



# 2 If vn is light during inflation, such a displacement is generically expected to arise from the quantum 
fluctuations of the vr. The fluctuations will become isocurvature fluctuations in the CDM. It is also 



possible to generate large amount of non-Gaussianity in the CDM isocurvature perturbations [42 1 
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abundance of the observed dark matter. Note that the Qp R is determined solely by the 
reheating temperature and the initial position of vr in this case, and is independent of the 
sneutrino mass. 



V. DISCUSSION AND CONCLUSIONS 

We have assumed that the i?-parity is dominantly violated by the bilinear term. If it 
is violated mainly by the trilinear terms, LLE, QDL, or UDD, the predicted cosmic-ray 
spectra will be different. If the i?-parity is broken by the LLE operator, the vr will decay 
into charged leptons as well as neutrinos through the left-right mixing, which will result in 
a sharp peak in the predicted positron fraction at the solar system. The continuum gamma 
rays are also produced, while virtually no antiprotons are produced. This choice therefore 
provides an ideal way to avoid the observational constraint on the antiproton flux. In the 
case of QDL or UDD, the hadronic branching ratios are larger than the case of LLE, and 
therefore we expect a relatively large contribution to the antiproton flux. 

With the tiny i?-parity violations in our scenario, the NLSP is generally long-lived in 



collider experiments and decays outside the detector [43|, 144 1. If NLSP is the neutralino, it 
is observed as a missing energy, and the collider phenomenology is the same as the i?-parity 
conserved case. If stau is the NLSP, we expect to observe its track inside the collider, which 



is similar to the case of the decaying gravitino LSP js|. However, there is one possibility 
that the vr is only a fraction of the total dark matter, and the rest is explained by some 
other stable particles, such as a QCD axion. If the fraction of vr in dark matter r is very 
small, the i?-parity violating coupling should be enhanced as oc r~ 1//2 , for a fixed positron 
flux. Since the lifetime must be longer than the present age of the universe, r can be as 
small as (9(10 -10 ), and the corresponding i?-parity violation can be enhanced by a factor 
of O(10 5 ). With such a large enhancement in the i?-parity violation, we may observe some 
signatures of decaying NLSP at colliders. Note also that, even if the fraction is much smaller 
than unity, the features of the cosmic-ray spectra from decaying right-handed sneutrino will 
not change. 

We have focused on the case that all the right-handed sneutrinos have almost the same 
mass, and considered only one of them is the source of cosmic-ray. If the three right-handed 
sneutrinos exist with an equal amount in our universe, they will all contribute to the cosmic- 



1(3 



ray. The contribution from one of the three right-handed sneutrinos may dominate over those 
from the other two, depending on the values of \ii and the neutrino mass spectrum. 

Even there exist mass differences among the three right-handed sneutrinos, unless the 
decay channel for the heavier one to a on-shell higgsino is kinematically allowed, the decay 
patterns of all the right-handed sneutrinos are similar to each other. Namely, the heavier vr 
also decays via i?-parity violating interactions, since the i?-parity conserving decay channel 
is suppressed as explained in the following paragraph. 

We should also emphasize here that the excess in the cosmic-ray positron can be also 
explained without introducing the i?-parity violation interactions. If the right-handed sneu- 
trinos are the three lightest supersymmetric particles, and if there is slight hierarchy in the 
masses, the heavier one will decay into a lighter one accompanied by two SM leptons through 
an off-shell higgsino. It is obvious that the two SM leptons will be the source for the cosmic 
positron and gamma-ray excesses. Also, the antiproton will not be produced. However, the 
lifetime tends to be much longer than the needed one to account for the positron excess, 
since the decay amplitude is suppressed by the neutrino Yukawa coupling squared. One 
solution is to add a small Majorana mass term for the right-handed neutrinos [jjj]. Then 
the see-saw mechanism will occur at very low energy, and the neutrino Yukawa coupling 
becomes larger, which can make the lifetime of the heavier one short enough to explain the 
positron excess. Another way to make the lifetime shorter is to tune the higgsino mass close 
to the heavier right-handed sneutrino mass. 

It is even possible to explain the positron excesses in a non-supersymmetric theory. One 
candidate is the sterile neutrino dark matter, whose decay produces charged leptons and 
neutrinos. The decay processes of a sterile neutrino lighter than 100 MeV were studied in a 
different context (see e.g. 45j). Since the sterile neutrino must be heavier than 100 GeV in 
our case, it will decay also into W and Z bosons, and also into wy with a suppressed rate. 
Thus, the energy spectra of the cosmic-ray particles from the sterile neutrino decay will look 
similar to those in the case of the decaying gravitino. The mixing angle between the sterile 
and active neutrinos must be extremely suppressed in order to realize the long lifetime of 
the sterile neutrino, which will make it difficult to produce the sterile neutrino through the 
mixing. It will be necessary to introduce a non-thermal production mechanism of the sterile 
neutrino. 

So far, only the normal hierarchy case in neutrino mass spectrum is considered, but we 



17 



expect that an essential feature of our results will persist in the other neutrino mass spectra. 
Rather, it may be possible to extract some information on the neutrino mass spectrum from 
the positron and gamma-ray fluxes by the PAMELA and FGST in operation. We leave 
those issues for future work. 

We propose a leptonic dark matter candidate whose decay may explain the excesses of 
observed cosmic rays. If the dark matter has a lepton number, it will not be surprising 
that we see some excesses only in the positron and the gamma-ray fluxes while virtually 
no primary antiprotons are observed. To illustrate the idea, in this paper we have studied 
a scenario that the neutrinos are Dirac type, and the right-handed sneutrino accounts for 
the observed dark matter and decays into the SM particles through the i?-parity violating 
interactions. The charged leptons as well as neutrinos are directly produced during the 
decay, leading to a sharp rise in the positron fraction. The decay products will also generate 
continuum gamma rays, which may account for the EGRET excess. Interestingly, with 
an appropriate set of the parameters, the antiproton production can be suppressed in our 
scenario, which is consistent with an observation that most of the observed antiprotons are 
secondaries. Furthermore, we expect that the energy spectra for those cosmic-ray particles 
may provide a new probe into the neutrino mass spectrum. Those features on the predicted 
cosmic-ray spectra in our study will be checked by the PAMELA and FGST satellites. 

Note added: After our paper was posted, the PAMELA group reported a steep rise in the 



positron fraction 46(, which can be explained by our model with a 300 GeV right-handed 



sneutrino dark matter, as shown in Fig. [3l 
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